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Abstract 

The asymmetry in the angular distribution of np — > d-K^ due to Charge Symmetry Breaking is 
calculated using Heavy Baryon Chiral Perturbation Theory. Recent developments in power count- 
ing have proven successful in describing total cross sections, and we apply them to the asymmetry 
calculation. Reducibility in one of the leading order diagrams is examined. We compare the up- 
dated theory with experimental results for a range of physically reasonable parameters and find 
over-prediction for the entire range. 
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I. INTRODUCTION 



The nuclear reaction A^A^ NNtt has been studied for a long time in many ways. 
Theoretical understanding of the reaction is currently pursued using Heavy Baryon Chiral 
Perturbation Theory (HB^PT) This effective theory of the strong nuclear 

force treats hadrons and mesons as the fundamental degrees of freedom. The Lagrangian 
contains an infinite number of terms which decrease in importance according to an expansion 
in a (somewhat) small parameter. Such an understanding is important because it is model 
independent, provided convergence is achieved. Furthermore, once theorists arrive at an 
ordering scheme for the expansion, error estimates will become more reliable and calculations 
of different and more exotic reactions possible. 



One particularly interesting observa 
of Charge Symmetry Breaking (CSB) 



ale that HB^PT can help determine is the magnitude 



m 

imate 

invariance of hadronic systems under an isospin rotation of vr about the y-axis. This sym- 
metry is broken by the mass difference between up and down quarks and by electromagnetic 
effects jsj. We will discuss the corresponding symmetry breaking terms in the HByPT 
Lagrangian. These terms are especially important for the interactions of neutral pions 9|. 

Relations involving the light quark masses can be obtained from the SU(3) chiral symme- 
try on which HB^PT is based, but experiments which directly measure — m„ are difficult 
because of the lack of a neutral pion beam. One experiment which overcomes this difficulty 
and also minimizes electromagnetic effects is np dn^. The angular distribution of this 
reaction is symmetric about 90° in the center of mass when charge symmetry is respected. 
Ref. [lo| recently observed that this distribution is asymmetric. 

This report advances previous work in several ways. The authors of Ref. llj] showed 



that a vertex which was thought to be higher order in fact contributes at leading order. 
This lead to a much improved understanding of the total cross section of pp dn^ . We 
extend this calculation off threshold for neutral pion production and make a comparison 
with the corresponding data. We also investigate a subtlety regarding reducibility in one of 
the leading order diagrams and comment on the resulting corrections. Finally, we calculate 
the asymmetry of np dir^, bringing up to date the calculations of Ref. 12 1 and 1^ by 
including the new effects discussed in this paragraph. 

In Sec. [TTl we discuss the kinematics and selection rules of the np —>■ drc^ process. Then, 
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in Sec. IIIII we develop the formalism necessary for calculation of the cross section. Here we 
detail the procedure for handling strongly interacting initial and final states. In Sec. IIVI 
we present the leading order diagrams with vertices from the isospin conserving part of the 
Lagrangian (Appendix . In this section we also provide a review of the power counting 



developed by Ref. and discuss its impact on neutral pion production. Next, in Sec. El 
we present diagrams with vertices from the CSB part of the Lagrangian which contribute to 
the reaction. Our results are given in Sec. |VT] and discussed in Sec. I VI II Also in Sec. IVIH 
we compare our work with another recent calculation of the asymmetry by Ref. jl^ . 

II. KINEMATICS AND SELECTION RULES 

At threshold, the reaction A^A^ — > NNtt produces a pion from the kinetic energy of the 
incoming NN pair. Let p = \{pi — P2) be the relative momentum of the pair. In the center 
of mass frame for this reaction, the total momentum P = pi + p2 = and thus pi = p and 
P2 = —p. In the non-relativistic limit, the total initial energy is Ei = 2Mn +p'^/Mn where 
Mtv is the nucleon mass (the average mass is used for np reactions). The energy of the final 
state at threshold is Ef = 2Mn — Ei, + where Et = 2.224 MeV is the deuteron binding 
energy and is the mass of the appropriate pion. Neglecting the binding energy yields a 



simple scale for external momenta in these reactions, p = \/r(^J<A^ = 356 MeV. 

There are two frames to keep in mind: the center of mass frame (C), and the lab frame 
where the proton is at rest (E for experiment). The experiment of interest [lO] is np dir^ 
performed in the E frame. In this frame the invariant is expressed as se = {Mp + M„)^ + 
2MpTL where Tl is the kinetic energy defined by En = Mn + T^. The experiment was 
performed at = 279.5 MeV, slightly above the threshold value of = 275.1 MeV. To 
simplify the formalism, we use the C frame to do the calculation. In terms of the pion 



momentum, q, and the deuteron mass. Ma, we have y/sc = a/ + g ^ + \/m1 + g^. It 
is convenient to define a dimensionless parameter rj to describe how far off threshold the 
reaction is, rj = |g l/m^r. Equating the invariants we find rjc = 0.169 (gc = 22.86 MeV) at 
the experimental energy. 

We use the hybrid approach developed by Weinberg [l5] to calculate the cross section, 
due to the fact that the initial and final states are strongly interacting. The "operator" is 
calculated as the sum of two-particle-irreducible diagrams involving four nucleon lines (two 
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TABLE I: Channels for the np wavefunction in np 





Strong 


CSB 




'Pi 


'Pi 




'So, 'D2 


''Si, 'Di, 



incoming and two outgoing) and one pion line. Then, a phenomenological potential is used 
to calculate the NN scattering wavefuctions and the deuteron bound state wavefunction. 
Finally, the operator is convolved with the wavefunctions to obtain the matrix element. 
This method is described in detail in Sec. Illlt but first we present the selection rules for 
np — s> dn^. 

The deuteron has the following quantum numbers: spin S" = 1, orbital angular momentum 
L = 0, 2 (parity even), total angular momentum J = 1, and isospin T = 0. If the pion 
{J'^ = 0~) is produced in the s-wave, then conservation of total angular momentum gives 
Jj = 1. Furthermore, since the pion is parity odd, s-wave pions must be produced from a 
parity odd np wavefunction, Lj = 1. Likewise, if the pion comes in the p-wave, the initial 
parity is even and Jj = 0,1,2. 

To completely pin down the quantum numbers of the initial state, we turn to isospin. 
If isospin is conserved in the final state, then Tf = 1 because the deuteron is an isoscalar. 
"Strong" diagrams have Tj = 1 with the isospin part of the initial wavefunction symmetric 
under exchange of nucleons. Overall wavefunction antisymmetry then requires Si = 1 for 
In = and Si = for 1^ = 1. Thus the initial np pair must be in the 'Pi channel for 1^ = 0, 
while the available channels become ^5*0 and for 1^ = 1. CSB operators transform as 
vectors under isospin, thus the initial neutron-proton state must have Tj = 0. The same 
arguments as for the strong operators then show that s-wave pions are produced from ^Pi 
np pairs while p-wave pions are produced from the coupled channels 'Si and 'Di in addition 
to 'D2. These conclusions are summarized in Table [H 

The observable of interest in the experiment is the forward/backward asymmetry in the 
differential cross-section given by 

A non-zero asymmetry will only be observed when initial states of opposite parity interfere. 



However, the interference can only occur for states with the same spin since the spin z- 
components get summed over. Thus for calculating the asymmetry, we are concerned with 
two terms: (s-wave strong) ■ (p- wave CSB) and (p-wave strong) ■ (s- wave CSB). 

III. CROSS SECTION 

A. Cross Section Formalism 

We will now derive an expression for the cross section in the C frame where the pion 
momentum is q. First we must define expressions for the strongly interacting initial and 
final states. To form an interacting NN state with total momentum P = pi + p^, we use a 
superposition of free particle states 



where spin and isospin have been ignored for now. The wavefunction il){p) is obtained by 
solving the Schrodinger equation with the appropriate NN potential specified below. In 



where pi_2 = (£'1,2? ±p), 1 = (^^g,^), K = {E^.—q). Note that we are treating the initial 
state as two separate particles, but the deuteron as a single particle. As mentioned in Sec. 



is calculated in momentum space and is a function of the external momenta p, fc, and q. Also 
note that the wavefunctions will include spin and isospin kets on which the operator acts. 

We perform the calculation in position space by inserting Fourier representations of both 
wavefunctions 






nil the sum of the two-particle-irreducible diagrams with external momenta q,pi,2, ^1,2 (the 
"operator") is denoted Ai, and is convolved with the external wavefunctions. The operator 




(4) 



A Fourier representation of the operator with respect to / = A; — p is also inserted 




(5) 
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As described in Appendix [Cl I is the momentum that appears in pion production reactions: 
Now we can rewrite Eq. ([3]) 



(27r)3 (27r)3 

d\rdir)M{f^ iPnpir). (6) 

With these choices of Fourier representations, the momentum integrals evaluate to delta 
functions which allow evaluation of the spatial integrals. The result is an integral over a 
single spatial variable. 

Next, we express the invariant S-matrix element in terms of Ai 

{'K\q)d{K) I S I p(pi),n(p2)> = l-i{27^f6\q + K-p^ ~P2)M{puP2 - q,K). (7) 

In the center of mass frame, the differential cross section is 

d^q d^K 



da = ^i E \Mn2.Y5\q + K-p,-p,)- 



(27r)32cj„ (27r)32Erf' 

m^, mi ,7712 y \ / 



(8) 



where we have averaged over the four np spin states and summed over the three spin states 
of the deuteron. The vector part of the delta function tells us g* = —K and the energy part 
tells us El + E2 = Uq + Ed- This allows us to perform all but the dflx integral, 

— = ^Li y \M\' (9) 

' ' m£j,mi,m2 

What remains is to derive expressions for the wavefunctions and the operator appearing in 
Eq. (p. 



B. Initial and Final States 

In the absence of interactions the wavefunction ipnpi^) is determined by performing a 
partial wave expansion on an antisymmetrized wavefunction of a free proton and a free 
neutron with relative momentum p. First we consider the strong operators where the np 
pair is in an isospin-1 state, 

^np) = PT=i^(e'^^1mi,m2)®|T,,i,T,,2)-e-*^^1m2,mi)®|T,,2,T,,i)) 

= (e^^" Vi, ^2) - e-*P" V2, mi)) ® |T = 1, = 0) , (10) 



r 



where Vt=i is the isospin projector. The bra (r | indicates that we are choosing a basis for 
space, but not for spin or isospin. Imphcit in the notation is the dependence of (r | ipnp) 
on the momentum p, the spin z-components of the two nucleons, mj, and the isospin z- 
components of the two nucleons, j, with the requirement that i + Tj. 2 = 0. 

The exponentials are now expanded and the presence of the strong interaction is added by 
changing the spherical Bessel functions, jiipr) e^^^ui^jir) /pr . The ul,j functions and the 
5l phase shifts are obtained by solviiig the Schrodinger equation with a phenomenological 
NN potential (we use Argonne V18 [l6|]). Finally, the spherical harmonics are combined 
with the spin kets to form states with definite total angular momentum. The notation for 
these states is \{SL)J,mj) ® \ T,Tz). For the allowed quantum numbers, we find 

(r I ^npC'^'Lj,T= 1)) = 47r(0^e^^-^^^ (1/2 m^, 1/2 \ S m,) 
X J](^m„Lmi-m, I J m,)Y'^^:„^(p)(f I (5L) J, mi) ®|1,0), (11) 

mi 

where = mi + m2 and the second Clebsch Gordan coefficient allows us to make the sum 
over rrii = mi +ms rather than m/. For the CSB operators, we have T = np wavefunctions 
and find 

(r I iJnpe'+'Lj,T = 0)> = ±47r(2)V^^^^^^ (1/2 mi, 1/2 ms | S m^) 
X ^ ( J I 5 m„ L m, - m,) Y^;_,^^ (p) (f | {SL)J, m,) ® |0, 0) , (12) 

mi 

where the ± refers to T^, 1 = ±1/2. Similar analysis gives the final-state wavefunction of the 
deuteron as a function of its polarization, rrif, 

(Mmf) I f ) = (0, 0| ® ^ ((10)1, mj\f) + ^ ((12)1, | f) . (13) 



IV. STRONG CONTRIBUTION 



A. Diagrammatic Expansion 



Before we calculate the effects of charge symmetry, we need to discuss the power counting 
scheme which leads to a calculation of the total cross section in agreement with experiment. 
HB^PT orders contributions in powers of the external momenta divided by the chiral sym- 



metry breaking scale, which is ~ M, 



N 



ll| . In this inelastic reaction, both q and p appear as 



(a) 



ib) 



FIG. 1: Leading order contributions to np — > di:^. Solid lines represent nucleons, the double solid 
line represents a A, and dashed lines represent pions. 

external momenta and we need to keep track of both in the power counting. We define the 
expansion parameter x = p/M^ = ^Jm-^/MN = 0.40. 

The Lagrangian is given in Appendix |Al The index of a "type vertex is given by 



T^i = di + ^ - 2 



(14) 



where di is the sum of the number of derivatives, m^'s, and 5's (the AN mass difference), and 
fi is the number of fermion fields. In standard power counting at tree level, the sum of the 
z/j for each vertex in a diagram indicates the power of x which that diagram contributes. 
This rule, however, will require modification due to the relatively large value of p. 

There are three two-particle irreducible diagrams which can be drawn using the vertices 
from C^^\ They will be referred to as the impulse (Fig. [T^), rescattering (Fig. [lb), and 
Delta (Fig. ^) diagrams. There is some ambiguity as to what is meant by Fig. [1^. We will 
discuss this in detail below. First let us make a few observations about these three diagrams. 
We have stated that the initial relative momentum is large. The final state nucleons (the 
deuteron) have a comparably small relative momentum. This "momentum mismatch" is 
provided for in the rescattering diagram, but not in the impulse diagram. For this reason, 
the rescattering diagram dominates the total cross section. This diagram is strongest in the 
channel with s-wave pions, and the ttttNN vertex is typically referred to as the Weinberg- 
Tomozawa (WT) vertex. Next, although the Delta diagram provides the momentum transfer 
required, the Delta resonance is at 1232 MeV and the vrA^ energy is ~ 1080 MeV, so the 
Delta diagram is also somewhat suppressed for our situation of interest. Finally, we note 
that both the impulse and the Delta diagrams are strongest for the channels in which the 
pion is in a p-wave. 

We now return to the issue of the impulse diagram. On the one hand, we know that a 
single nucleon cannot emit a pion and remain on shell. But on the other hand, the diagrams 
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(a) 




ib) 



FIG. 2: Kinematically correct, but reducible impulse contributions. Solid lines represent nucleons 
and dashed lines represent pions. 



FIG. 3: "Hybrid" approach. Solid lines represent nucleons, dashed lines represent pions, dotted 
lines represent pions with u; = 0, and filled ovals represent NN strong interactions. 

in Fig. [2] which remedy this problem by including one pion exchange (OPE) appear to be 
two particle reducible. To resolve this issue we first note that near threshold the energy 
of the exchanged pion in each of these diagrams is approximately uj ~ 771^^/2. However, 
the diagram of Fig. [1^ is evaluated as diagrams of the form of Fig. [2l This is due to 
hybrid nature of the calculation; once the operator (traditionally, the irreducible diagram) 
is calculated, it is convolved with NN wavef unctions. One of the major terms of the strong 
interaction potential at low energy arises from static OPE (cj = 0). We ignore the effects of 
the rest of the wavefunction for the moment. The effects of static OPE are schematically 
shown in Fig. [3l 

To obtain the correct impulse contribution, we add up the contributions from Figs. [1^ 
and [2] and then subtract what is already included in the wavefunctions (the last two diagrams 
of Fig. [3]). This calculation is schematically shown in Fig. |H 

The OPE propagator is of Yukawa form. 



DJr) 



(15) 



where fi{uj) = ^/m^ — uP' ^ \f%l1mT^. Thus subtracting off the final two diagrams in Fig. H] 
amounts to making the replacement 



-\/3niir''/2 



(16) 



FIG. 4: Complete impulse contribution. Solid lines represent nucleons, dashed lines represent pions, 
dotted lines represent "wavefunction" pions with uj = 0, and the square represents the operator 
which is used for the full impulse approximation. 

in the exchanged pion propagator. The final four terms of Fig. H] comprise a correction to 
the impulse diagram, which we will refer to as a "wavefunction correction." We find that 
this correction is ~ 4% of the total impulse amplitude at the experimental energy, and we 
include it in our calculation. 



B. Power Counting 

Now we will look more closely at the size of these diagrams in the "p kinematics" , using 



the counting techniques developed in Ref. 
Lagrangian in Eq. (lAip 

Dn{p) 



17l |. The propagators are calculated from the 



Da{p) 



p° + ie 

i 



— (5 + ie 



— ^ 

+ [m^ — [p^)^) — te 

Both the ttNN and the tcNA vertices have a momentum dependence of |g |, the pion mo- 
mentum at that vertex. Note that this momentum is p in the OPE verticies. The WT vertex 
contains a factor of Ug^in + ojq^out- 

The external particles have the same momenta in each diagram. The produced pion has 
q = {ijOq,q) ~ (m7r,0), and the incoming nucleons have pi^2 = (-Ei,2,±p) ~ ^"^7r/2,±pj 
in the Heavy Baryon formalism in which the nucleon mass is subtracted off of the energy 
component. Consider the impulse diagram of Fig. [2ti. The final emission contributes g, the 
nucleon propagator 1/^,,^, and the OPE p- 1/p^ -p, so that the diagram is ~ q/rriT^ = rj. The 
rescattering diagram is ~ 3?t^72 ^ ^j^^ Delta diagram is ~ Finally, note that 

?7 ~ and 6 ~ 2mjr so that while the rescattering diagram is ~ x the impulse and Delta 
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(a) 




ib) 



FIG. 5: Recoil corrections. Solid lines represent nucleons, dashed lines represent pions, and the 
filled circle and square represent = 1 vertices. 

diagrams are numerically ~ x^- This ordering comes in agreement with the fact that the 
rescattering diagram is strongest for s-wave pions while the impulse and Delta diagrams are 
strongest for p-wave pions. 

It is well documented that these three diag rams alone do not correctly reproduce the 



experimental data for the reaction; see Ref. 18|] for a review of the theory of meson produc- 
;ion. Near threshold [t] ^ 0.05), the most recent experiment found aexp{np —>■ dn^) ^ 90 /ib 
lol . while for these first three diagrams, we find a ~ 55 /ib. The solution to this problem 
was discovered by Ref. llj, who noticed that the u = 1 "recoil" correction to the WT 
vertex, which is found in Eq. ( 1A2I) goes like (gi„ + (jout) ■ {Pin + Pout) / i^M^) where g is the 
pion momentum and p is the nucleon momentum. For Fig. [5|l, this vertex (a filled circle) 
~ p /{2Mn) = and thus this diagram is of order x, the same order as the i/ = 

rescattering diagram. Similarly, we find that the s-wave portion of the recoil correction to 
the impulse diagram (Fig. [5)d) is of order x- In this diagram, the filled square represents the 
sum analogous to Fig. H] for the recoil diagram. We find that the wavefunction corrections 
are more important (~ 20%) in this case. Finally, the s-wave portion of the Delta diagram's 
recoil correction is found to be higher order and is therefore ignored. 

The recoil corrections to the propagators have also been included in the calculation where 
applicable. For this reaction, the only such diagram is Fig. |2]d where the 3-momentum in 
the nucleon propagator is large (~ p). For that propagator, we use the corrected version. 



Dn{p) 



(18) 



Using this propagator rather than the z/ = version doubles the size of Fig. Nevertheless, 
this diagram (minus its cUtt = analog) is already very small. Thus the net effect of correcting 
the propagators is small for this reaction at this order. 

Including all these recoil corrections (especially Fig. [5^) brings the theoretical cross 
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0.05 



0.10 



0.15 



0.20 



FIG. 6: Cross section for np dir^ in terms of a = a/rj as a function of the pion center of mass 

□ 

momentum, r] = q/m^o. Circles with error bars display the experimental results of Ref. [19[]. The 
dashed line displays the results of including the diagrams in Fig. [1] and the solid line displays the 
results of also including the recoil terms discussed in the text. 



(a) 



(b) 



id) 



FIG. 7: Irreducible loops. Solid lines represent nucleons and dashed lines represent pions. 

rn 

section near the experimental results as shown by the solid curve of Fig. [6] |19i]. Due to 
the relative scatter of the data shown in Fig. [6] it is difficult to tell how well the theory is 
reproducing the experiment. Regardless, it is clear that theoretical improvement has been 
made. 

It should also be noted that the subtlety of reducibility and recoil corrections in this 
reaction resolves questions about NLO loop diagrams discovered by Ref. [2^. Namely, the 
sum of all the NLO irreducible loops in Fig. [7] is found to be proportional to p. This is a 
problem because such sensitivity of the operator to the NN wavefunction is not physical. 
Again, we consider including OPE in the operator, this time for the rescattering diagram. 
There are two resulting diagrams shown in Fig. [SI It was shown in Ref. ll|] that in these 



reducible loops, the recoil corrections to the nuc 



addition to the WT's recoil correction. Ref. 



eon propagators need to be included in 
ll\ then showed that (part of) the energy 



dependence of the WT vertex "cancels" one of the nucleon propagators leaving a reducible 
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FIG. 8: Reducible loops. Solid lines represent nucleons and dashed lines represent pions. 

diagram similar to Fig. [T^i. This diagram is equal in magnitude and opposite in sign to 
the aforementioned NLO sum, resolving the issue. The other term that remains from the 
original loop integral after this manipulation is still of reducible form but now has an on-shell 
WT vertex ~ 2m.^. This term would already appear upon convolution of the rescattering 
diagrams discussed above (including recoil corrections) with external wavef unctions, i.e. this 
term is truly reducible. This result can be stated another way: cancellation of the irreducible 
loops comes from a short range piece of including OPE in the rescattering operator. Thus 
a complete NLO calculation must include this OPE. However, its remaining reducible piece 
would only be included in the calculation if the recoil correction to the WT vertex is included 
in the rescattering operator. Therefore it is both consistent and necessary to include the 
WT recoil correction which was shown above to reproduce the total cross section. 



P-Wave Observables 



Another important test of the theory is how well it describes p-wave pions [2l(]. This 
is especially important for the asymmetry, which involves strong p-waves at leading order. 
The differential cross section can be expanded in Legendre polynomials, 

^ = ao + aiPi(cos(^)) + a2^2(cos(^)) + (19) 

where 6 is the angle between p and q. Note that the total cross section plotted in Fig. [6] is 
a = AiraQ/r]. As discussed in Appendix [Dl a2 receives contributions almost exclusively from 
p-wave pions. The ratio a2/«o is therefore used as a test for this part of the theory. We 
find that the diagrams of Fig. [1] along with their recoil corrections overestimate the data by 
approximately a factor of two. Upon closer inspection we find that the ^5*0 amplitude (which 
is known to be small) is relatively large. This amplitude is coming mainly from the Delta 
diagram, as can be seen in Appendix [D] where the values of the reduced matrix elements are 
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Q'2/Q'o 




FIG. 9: Legendre coefficients of the differential cross section for different values of the cutoff. Data 



is from an np cItt^ experiment (Ref. .19|, circles) and an pp — > dir^ experiment (Ref. 23], 
squares) in which the data have not be corrected for Coulomb effects. 

listed. To remedy the situation in the simplest way possible, we implement a cutoff for the 
Delta diagram 



p"^ + jj? p ^ + /i^ \P^ + 



(20, 



+ /i^ p"^ + A? J — p? 

One can show that doing this essentially softens the OPE potential for r < log(A//i)/A. 
Clearly this is not an acceptable long-term solution for an effective field theory, but the fact 
that such a procedure is necessary is interesting given that the reaction occurs at an energy 
~ 150 MeV below the Delta resonance. That p-wave pion production is highly sensitive to 
the strength of its contact term was discussed by Ref. 22 1. 



Figure [9] shows the effects of this cutoff on the ratio 0:2/0:05 where A = 10 GeV represents 
the original theory (such a large cutoff has no significant effect). Note that the amplitudes for 
np dn^ are related to those for pp — > dn^ (which are bigger by ^/2) when charge symmetry 
is respected. Thus the ratio plotted should have the same value for both reactions. Also 
note that for simplicity we used a single cutoff, and that if we had used one at each vertex 
of the OPE we would have found A — > y^A. By adjusting the cutoff to fit the data, we find 
A = 310 MeV. 

Another useful observable for testing p-wave pion production is the analyzing power. Ay, 
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FIG. 10: Analyzing power for different values of the cutoff. Data is from a pp ^ d-Tr"*" experiment 



(Ref. 



circles) in which the data have not be corrected for Coulomb effects. 



which is defined 

_ da^{e)-dai{e) 
^'^^^ = da^{e) + dam ^ ^ 

^^^'^^^^ " eJl^i ^ |^(m,, = ±l/2,^)r, (22) 

where rrii^y = ±1/2 refers to the fact that the beam is polarized perpendicular to the 
scattering plane. In the z-basis, these states are 

M = ±1/2) = M(m, = l/2)±iM(,n,^-l/2)^ 

v2 

As shown in Appendix [Dl Ay is proportional to the product of s-wave and p-wave amplitudes. 
Figure [TO] shows the effects of the cutoff on this observable. Again, charge symmetry implies 
that Ay should be the same for both neutral and charged pion production. We find the best 
agreement with the data for A = 560 MeV. Below, we will display our results using both 
the original theory and a cutoff taken to be the geometric mean of two fits, A = 417 MeV. 



V. CHARGE SYMMETRY BREAKING 



The fact that the up and down quarks have different mass is reflected in the Lagrangian 
by including terms which break chiral symmetry [24]. The leading such terms are given in 
Eq. (lA2p and have coupling constants Stun and (5m„ which are constrained by 

SrriN + SniN = M„ - Mp. (24) 
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FIG. 11: Leading CSB contribution. Solid lines represent nucleons, dashed lines represent pions, 
and crosses represent 1^ = 1 CSB vertices. 



r 



X 



(c) 



FIG. 12: = 2 CSB contributions. Solid lines represent nucleons, double solid lines represent A's, 
dashed lines represent pions, crosses represent u = 2 CSB vertices, and the boxed cross represents 
the full impulse CSB diagram including OPE. 

The Sitln term has its origins in the quark mass difference and its size is ~ (m^ — m„) = 
e{md + rriu) with e ~ 1/3. The formahsm of xPT teUs us that (m^ + m^) oc m^, and 
so dimensional analysis along with the natural QCD scale, Mat, yields drriN ~ em^/Misi. 
The term is of electromagnetic origins, but is of the same order as 5mi^. These CSB 
operators appear in the rescattering diagram depicted in Fig. [TT] where the CSB vertex 
is denoted with a cross. The size of this diagram is 5mN/p ~ em^/ {MnP) = ■ Note 
that although the full nucleon mass difference appears explicitly in the Lagrangian at this 
order, the corresponding operator {N'^t^N) does not change the parity and thus does not 
contribute to an asymmetry. In Sec. IVIIl we will discuss another source of CSB coming 
from a more detailed evaluation of the diagram of Fig. [T|d which was made by Ref. 14j ]. 

Another CSB term given in Eq. ( 1A3P involves one derivative and one {f3i ~ tm^/Mf^) 
and is thus a z/ = 2 vertex with momentum dependence |g |. This vertex appears in the dia- 
grams of Fig. [T2] whose sizes are /?ig/m^ ^ e^X^- Fig- fT2b . the boxed cross represents the 
sum analogous to Fig. H] for the CSB impulse diagram. Again, the wavefunction corrections 
are small (2%). 

As mentioned in Sec. [Tll contributions to the asymmetry come from interference terms 
between the p-wave part of the strong amplitude and the s-wave part of the CSB ampli- 
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FIG. 13: Interference terms for the asymmetry in np dir^. Solid lines represent nucleons, double 
solid lines represent A's, and dashed lines represent pions. The filled circle and square represent 
u = 1 strong vertices and crosses represent CSB vertices. The boxes represent full impulse diagrams 
in the sense of Fig. [H 

tude, and vice versa. The issue is somewhat complicated because, in contrast to threshold 
emission, each diagram can contribute in both the s-wave and the p-wave. However, contri- 
butions to the sub-leading parity (s-wave for the impulse and Delta diagrams and p-wave for 
the rescattering diagrams) are formally higher order. For example, the strong rescattering 
diagram for p-wave pions comes with a from the WT vertex, a 1/p^ from the pion prop- 
agator, and a q from the pion emission in the p-wave. Thus the diagram counts as ~ rjx^- 
The contributions to the asymmetry are depicted in Fig. [131 Fig- [T3b includes strong p- 
waves and CSB s- waves and has size [r]) x (ex^). Fig. [T3b includes strong s- waves and CSB 
p-waves and has size (x) x (er^x^). Fig. [T3b includes strong p- waves and CSB s- waves and 
has size (rjx^) x {^X^)- Thus we find that in these kinematics Fig. [T^ (~ erjx^) is the LO 
contribution, and Figures [T^.c (~ erjx^) both come in at NLO. Other interference terms 
involving these diagrams are higher order, ~ trfx^ o^' smaller. Finally we note that this 
work is not intended to be a complete NLO calculation as loops and higher order vertices 
may come in at NLO. 
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VI. CSB RESULTS 



In this section we discuss the CSB results of our calculation. The coupling of the WT 
vertex (and its recoil correction) is determined by chiral symmetry and we use = 91.9 
MeV. For the impulse vertices we use the values qa = 1.267 and = 2.1(7^. We use the 
following masses: Mn = {Mn + Mp) /2 = 938.919 MeV, M„ - Mp = 1.293 MeV, Ma = 1232 
MeV, and m^r = mT^o = 134.977 MeV. Consider now the CSB coupling constants. The 
Cottingham formula can be used to obtain 5mjv = —(0.76 ±0.30) MeV 25|. The constraint 

lat 



13| 



of Eq. ( !24|) then fixes (5mjv = 2.05 ± 0.30 MeV. The combination of these parameters t 
appears in the asymmetry is 5mN — Smiy/2 = 2.4 ± 0.3 MeV. It is also noted in Ref. 
that other models predict different values for drriN leading to 1.83 < < 2.83. 

Even less is known about (3i, the impulse CSB coupling. As a starting point, Ref. 
notes that this CSB operator can be viewed as arising from vr — 77 mixing. The result shown 
is that 



As discussed in the review 



Also, Ref. 



0.10 < < 0.51. (26) 



27| gives {tt^\H\7]) = -0.0039 GeV^ and we use = 547.51 MeV. These values 
result in —0.47 < c,, < —0.21. Thus it is at least plausible that the (3i term could originate 
naturally from rj — it mixing. We note that the rj' could also give such a term, but do not 
consider it here. Using Eqs. ( l25l) and ( l26l) . we obtain —3.2 x 10"^ < /9i < —1.4 x 10~^. 



Note that the value used in the original calculation of the asymmetry by Ref. [12| was 
(3i = —8.7 X 10^^, which we refer to as the "extreme value". However, recall from the above 

2 

discussion that the natural size of the p-wave CSB is (3i ~ —tjt ~ — 6 x 10~^. Thus even 
though its origins may not lie exclusively with the rj, the aforementioned "extreme" value 
for j3i is not extreme at all from the effective field theory's point of view. 

In Table [Tll we display our results which are obtained by using Eq. to eliminate 

SniN so that the diagram of Fig. [11] is proportional to 

Sm^-^ = — (27) 

Because the asymmetry is linear in the CSB amplitudes (and therefore the CSB parameters). 



18 



TABLE II: Asymmetry in np dn as a function of CSB parameters 5mi\i and "LO" and 
"NLO" represent the sums discussed in Fig. [T3l The term "z" arises from the influence of the term 
proportional to Mil^^£p^ 





Aff,i6mN = 2.05 MeV, 
f3i = -3.2 X 10~^) X 10^ 


Order 


Delta Cutoff 


X 


y 


z 


LO (no recoil) 


None 


33.7 





-14.5 


54.6 


LO (no recoil) 


A = 417 MeV 


27.6 





-11.9 


44.7 


NLO (no recoil) 


None 


37.6 


1.4 


-16.2 


56.4 


NLO (no recoil) 


A = 417 MeV 


32.5 


1.8 


-14.0 


47.0 


LO 


None 


25.0 





-10.8 


40.4 


LO 


A = 417 MeV 


18.9 





-8.2 


30.7 


NLO 


None 


28.1 


1.4 


-12.1 


41.0 


NLO 


A = 417 MeV 


22.1 


1.6 


-9.5 


30.7 


NLO 


A = 310 MeV 


20.1 


1.7 


-8.7 


27.3 


NLO 


A = 560 MeV 


24.0 


1.6 


-10.4 


33.9 



we are able to present our results as a set of coefficients, {x, y, z} defined by 

Af,xW' = x- (^^^ + y ■ (A X 10^) + z. (28) 

The primary advance made in this work over the the previous calculation of Ref. pJi] 
(including the rescattering and impulse recoil corrections) is shown in moving from the top 
four rows to the next four rows. At LO this simply increases ao, but at NLO it affects both 
the numerator and the denominator of the asymmetry. 

The experiment of Ref. lOj found Afi, = [17.2 ± 9.7] x 10^^. The first calculation of 
the asymmetry used a A^A coupled channel formalism and included the CSB impulse vertex 
as well as other, smaller effects arising directly from the neutron-proton mass difference 
[l^ . This study reported Afb = —28 x 10""^. The second calculation included only the 
CSB rescattering vertex, and found Afb = 60 x 10~^ [l3^. Both these calculations were 



preformed before the work of Ref. Uj which brought the total cross section into agreement 



with experiment. Our work brings the asymmetry calculation up to date. 
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As shown in Table |TT1 the value of A fb is over-predicted for the set of physically reasonable 
parameters used in the last column. Nevertheless, for the most "extreme" set of parameters 
discussed above {Sm^ = 1.83 MeV, (3i = —8.7 x 10^'^) the cutoff NLO calculation yields 
Afh = 16.9 X 10~^. The effects of using different values for the cutoff can be seen in the last 
two rows of the table. 



VII. DISCUSSION 

We have mentioned that it is difficult to tell how well the current theory is reproducing 
the total cross section. We have also seen that there is reason for concern regarding the 
theoretical description of p-wave pions, which comprise the entire strong contribution to the 
LO asymmetry. Because the total cross section is dominated by the rescattering diagram, 
small changes to the p-wave amplitudes are able to significantly modify the asymmetry while 
only slightly changing the total cross section. As a temporary solution, we implemented a 
cutoff in the Delta diagram and thereby achieved acceptable agreement with the p-wave 
data. Another solution to this problem is to use a coupled-channel A^A potential for the 

n 

initial state. This approach was taken by Ref. [IJ] who were able to achieve good fits to these 
data without a cutoff, since the OPE of the Delta diagram is then part of the wavefunction. 

There are other significant differences between Ref. |l^ and this work. First of all, they 
discovered an additional CSB contribution to the asymmetry which is a consequence of the 
time derivative in the WT vertex. The effect of this contribution is equivalent to a change 
in the CSB rescattering diagram. 



(29) 



2 2 2 

Thus in order to update our calculation, we drop the fifth column of Table [TTl (the "z" 
column). Secondly, they used experimental data (from pionic deuterium) to determine the 
Legendre coefficient, ao = 1.93 fih. This is significantly larger than the theoretical value we 
use, ao = 1.49 /ib at LO (ao = 1.28 /ib for A = 417 MeV), and leads to a smaller value 
for the for the asymmetry. Note that experiments for neutral [19] and charged 23|| pion 
production found ao = 1-39 /ib and ao = 1.64 fih (Coulomb corrected) respectively. Finally, 
they do not include the NLO CSB diagrams. 

For the sake of comparison, we used our code to calculate the LO asymmetry with 
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A = 417 MeV, without the "z" column, using ao = 1.93, and using their quoted values for 
gA and For these choices we obtain Afb = 14.0^^ x 10"'', which is to be compared 
with their result of Afb = 11. 5^^- x 10"^. They use the Cottingham sum rule to obtain 
6mN = 2.0 MeV and thus Aft = 23 x 10~^, which only overestimates the data by 0.6a. 
Finally, we display our best result without the "z" column, 

A^, = 22.1^ + 1.6(Axl03). (30) 

For the set of parameters used in Table HTl Afb = 40.2 x 10"'^ which is an overestimation of 
the data by 2.4 a. 

Several issues remain to be understood theoretically. Firstly, it appears that a large 
contact term will be required to suppress the ^5*0 channel in the strong amplitude if one uses 
a purely NN initial state. The interesting physics observation here is that the Delta part of 
the NN wavefunction seems to be much more active than it is NN scattering. Secondly, the 
experimental cross section used by Ref. 13] is not well predicted by theory at NLO, and 
this fact plays a large role in the overprediction of our calculation. This situation becomes 
even worse when a cutoff is used to decrease the p-wave amplitudes. Finally, the existence 
of multiple mass scales greatly complicates the power counting for this reaction, and it is 
clear that a converging expansion cannot yet be definitively claimed. For these reasons, 
we conclude that further calculations are necessary. In particular, one should extend the 
calculation to next order while examining both the power counting of recoil terms and the 
reducibility of loops. 
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APPENDIX A: LAGRANGE DENSITIES 



The 1/ = lagrangian of HB^PT (with isovectors in bold font) with the Delta included 
as an explicit degree of freedom is j28| 

- -^N^ (r ■ (tt X -k)) N + I^A^t ('r ■ a ■ Vtt) iV 
+ {ido -6)A + ^\nUt-S- Vtt) a + h.c. 



+ ..., 



(Al) 



where T and S are the 3/2 1/2 isospin and spin transition operators, and r and a are 
the pauli matrices acting on the isospin and spin of a single nucleon. The "+..." indicates 
that only the terms used in this calculation are shown. 

The V = 1 lagrangian includes propagator corrections, recoil terms, and the leading 
s-wave CSB operator 



£(1) 



2m N 
1 



2m 



N 



2m 



N 



+ 



1 2hA 
2mN 2/^ 



-X^iN'^T ■ (tt X Vtt) ■ VA^ - ^iNW ■ ■ka ■ VN + h.c. 

4f2 V J 2/, 



iN^T'-kS-VA + h.c. 



AT +5^^ Art 



T-3 + 



T ■ TT 



Ar + ....(A2) 



Although there are a host of z/ = 2 terms, we list the just the CSB term relevant for this 
calculation 

Pi .ri^ .r : (A3) 



/:(2) = i^jVta.Vvr3Ar+.... 



APPENDIX B: DEFINING REDUCED MATRIX ELEMENTS 



One can show that for s-wave pions, the production operator is always either a scalar 
or a rank-two tensor while for p-wave pions, it is always a rank-one tensor. This guides 
the following definition of the reduced matrix elements. Note that the Clebsch-Gordan 
coefficients which will be summed over as well as the Spherical Harmonics describing the 
angular distribution of the differential cross section are "pulled out." First we define the 
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strong reduced matrix elements, 

(/ I MtlU I .) = (ij^. + (""/.^^ll '"/) ^,) 

X (1 rris, 1 m/ - m, I 1 m/) (p) (Bl) 

(/ I I ^) = ^i?e^^ocgo5^„oyo°(p) + (2^/-l^|l^/) ce-^.cgoy^;(p), (B2) 

where = mi + m2, ego = (1/2 mi, 1/2 m2 | 0), and cg^ = (1/2 mi, 1/2 m2 | 1 m^). Aq 
and are the s-wave reduced matrix elements, and B and C are the p-wave reduced matrix 
elements. To clarify the notation consider A2, for example. 

A^^j irr^ (!i^((10)l II +^((12)1 ||) All,,., (4..^ || (11)1>) , (B3) 

where the subscript J = 2 on the Ai indicates that we are using the portion of the operator 
that is a rank-two tensor in the space of total angular momentum. Note also that we have 
used the following general definition of a reduced matrix element, 

{{S'L')J'm'j I 7j I {SL)Jmj) = ^'^'^'^ ^j, ^^'^ -{{S'L')J' \\T^ \\ {SL)J). (B4) 

Similarly for the CSB reduced matrix elements, 

(/ I Mt^ I .) = + e*cg„r-tf) (B5) 

(/ I MfX, I i) = 'J^' (i?„e-'" +B,e"») cg,t,.„.y„°(p) 

^ (lmi,i|2!!V)(c„e-'-+C,e-'.)cg, 

X (1 m„ 2 m/ - m^ I 1 m/) Yl*_„^^ {p) 
+ (2 m/, 1^1 ^^f) De^6,^g,{lms,2mf-ms \ 2 m^) (B6) 

where Ao and A2 are the s-wave reduced matrix elements and B, C, and D are the p-wave 
reduced matrix elements. Also note that the strong phase shifts have been denoted Sl for 
each of the three initial channels, the CSB phase shifts are denoted 6l for the ^Pi and 
channels, and the coupled channel phase shifts are and 6 p. Since the ^5*1 and ^Di channels 
are coupled, B and C are split into a and (3 parts which have different phase shifts in the 
presence of initial state interactions. 
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ki = (fcf , k - q/2) ki = (fc", ~k- q/2) 



3^. 



g' = K/2,g') 



Pi = (p'\p) P2 = (p", -P) 

FIG. 14: Strong rescattering diagram. Solid lines represent nucleons, dashed lines represent pions, 
and the pions' isospin z-components are 3 and a. 



Finally, for comparison pur poses we include a translation between our reduced matrix 



elements and those of Ref. 



2l| 



Co 
Ci = -i 

C2 



1 



1 



'Qtt 



-.Be 



iSi 



Arc 



(B7) 



where the "C's" are pp — rfvr"^ amplitudes and isospin symmetry has been used to determine 
the translations. 



APPENDIX C: DIAGRAM TECHNIQUE 

To establish the diagram technique consider Fig. [THin light of Eq. ([5]). We write down 
the amplitude using momentum space Feynman rules, Fourier transform, and then convolve 
with the initial and final state wavefunctions. The left line is taken to be "nucleon 1" and 
we make the approximation that the exchanged pion carries half of the produced pion's 
energy, q' = {ujq/2,q'). Momentum conservation gives q' = k — p + q/2. According to 
Eq. ( lAB the WT vertex contributes 1/(4/^) ri,6 easb (u;q/2 + Uq) while the nNN vertex 
contributes qa/ i'^fn) T2,a ^2 ■ {—q')- The momentum space propagator is —i/^q'"^ + yU^) where 
/i^ = — C(jg/4 is the effective mass of the rescattered pion. Next, as discussed in Sec. 
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nil At we Fourier transform with respect to I = k — p, 



e — — — = e " ' / - — ^ e 



,2 



(27r)3 g'2 + /i2 J (27r)3 q'^ + fi 

Anr 

= 'J^e-"^-''/n{r)a2-r, (CI) 

where h{r) = (1 + l//ir)e~^'"/r. 

The deuteron has isospin and the np wavefuction includes a T = 1, = isospinor 
|1, 0), and thus 

(0, I tea3bri,br2,a I 1, 0) = -2. (C2) 
At this point, we have (defining M' = Mj^J^E^ 2^2 2^^) 



0, I I 1, ) = — ^ Kr) e-"^'"'^ <?2 ■ r. (C3) 



To calculate the diagram with the WT vertex on nucleon 2 we consider how each part of 
the left side of Eq. flC3l) transforms under 1^2. Since the strong part of the Lagrangian 
is invariant under isospin, M. is invariant. The initial isospin ket, |1,0) is invariant as well, 
but |0,0) — |0,0). Also note that r — > — r. Thus, 

0, I M'^ I 1, O) = -^m^ Kr) e--^/^ a, ■ r. (C4) 

Defining 5* = (ai + (T2)/2, S = (cti — (T2)/2, and 

£^ = exp(zg ■ r/2) + exp(— ■ r/2) 

O = exp(zg ■ r/2) — exp(— ■ r/2), (C5) 
we have the complete rescattering operator, 

^0, I M'^s I 1, O) = -I ^Rs h{r) 5 ■ f + O S ■ f) , (C6) 



where 

To proceed, we preform a partial wave expansion on £ and O and just keep the leading 
term. Note that we use the coordinate system defined by g = i. Then we calculate the 
spin-angle matrix elements of the rank zero £ S ■ r and the rank one (9 S ■ f operators. The 
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final expression for {f \ Ai \ i) simplifies to a radial integral which is computed numerically. 
Note that the first term in Eq. flC6l) corresponds to s-wave pions because S carries L = 
while f carries L = 1, and thus the operator will change the parity. Likewise, the second 
term corresponds to p-wave pions. In terms of the reduced matrix elements of Eqs. fIBip 
and (lB2p . we have 



Af = 



= ^/2E,2E22Ed8n^RsV3Ko 



= ^2^1 2^2 2Erf 877 7^5 76/^2, (CJ 

where the integrals are defined 



Kq = I drr^ { - 2 



A/2r / V 2 / pr 



).(f)M.) 



-\/2r / V 2 / pr 

A,. / ( !^ - 2^) (I) M,.)-;^. (C9) 



APPENDIX D: OBSERVABLES 

One experimental observable is the analyzing power, Ay, defined in Eq. (I?T|) . In the 
strong sector, we find (neglecting A2 which is numerically small) 

^ _ V3 cos(0) sin(g)Ao {V2B sm{6o - 61) + C sm{62 - ^1)) 
^ ~ 3^2 + 52 + C2 + (C2 - 2v^5Ccos(52 - 5o)) ^2(cos^) ' 

where the angular dependence is that of the nucleon relative momentum, p, with respect to 
the pion momentum, q = z. To compare with experimental results, which use p = z and 
(pTT = 0, we need to set (j)N = tt and 6^ = 9.„. 

To calculate the differential cross section as well as the asymmetry, we need to square the 
sum of all the matrix elements, sum over rrif and average over mi and m2- First we define 

= Mo + Ml Pi (cos 0) + M2P2(cos^) + M3P3(cos^), (D2) 
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so that 









64 TT^ S 


p\ 



AtiMq 



A 



Ml - IM3 



fb 



2Mn 



(D3) 
(D4) 



The results for the required quantities are 



Mn 



Ml 



24 TT 



3{aI + aI) + -{aI + aI) + ^b' + bI + bI + c' + cI + cI + d' 

+ 2 (B^B^ + CaCp) cos {6a - 6f,)] 
B ^0 - 2y^^2 j cos(5i - 5o) - V2C (ao - ^^^2 j cos(5i - 62) 



(D5) 



Ba + "^C*" ) cos((5c, - 61) + {a^ p) 



SVIOtt 



- A / ^ (^Ao - ^4j^^2 ) ^ C0S((52 - 5i) 



{A0A2 - 2A0A2) - -^{Al - 2aI) 
zv iU 



(D6) 





(25Ccos((52 - (5o) - - BpCp - (B^Cp + BpC^) cos(5« - 



^ (^Ba + ^C*") ^ C0S(52 - 5a) -(«->/?) 



CA2 cos(5i - ^2) - -^A2D cos((52 - Si) 
v3 



(D7) 

(D8) 



AOtt V 5 

Disregarding the small M3 term, the asymmetry is proportional to Eq. (lD6p . The physical 
content of Eq. flD6P 's first line is the interference of strong p-wave pions with CSB s-wave 
pions and the second and third lines are strong s-wave and CSB p-wave. 

Table HIT] shows the strong reduced matrix elements and Table HVl the CSB reduced matrix 
elements. The CSB rescattering numbers were calculated including the new contributions 
discovered by Ref. [l4|. The np phase shifts (in radians) which appear in the cross section 
according to Eqs. ( 1D5IID8I) are given in Tables |V] and IVIi 
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TABLE III: Strong reduced matrix elements 



Diagram 


^0 


A2 


B 


C 


Impulse (w/ wfn corr) 








-6.59 


32.85 


Impulse Recoil 


-5.62 


-0.21 


1.17 


-8.56 


RS (w/ Recoil) 


81.12 





0.54 


1.66 


Delta (no cutoff) 








-33.94 


37.96 


Delta (A = 417 MeV) 








-10.48 


21.60 



TABLE IV: CSB reduced matrix elements 



Diagram 


^0 






Bp 


Cct 


Cp 


D 


Impulse (X;^;-) 








12.23 


29.72 


-7.80 


-15.05 


-28.30 


RS(xm^) 


-28.83 





-1.37 


1.79 


1.48 


1.89 


-4.92 


Delta (x-i) 








12.60 


-7.71 


-8.47 


-2.86 


22.37 
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